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摘要 DA 白矮星 的 光谱 在 光学 波段 主要 由 巴尔 默 线 主导 。 由 于 谱 线 比 较 宽 ， 而 且 谱 线 轮廓 不 对 称 ， 传 统 
的 线 心 确定 视 向 速度 方法 非常 困难 。 本 文 介 绍 了 一 种 基于 利用 白矮星 的 有 效 温 度 和 表面 重力 加 速度 来 选择 
理论 模板 通过 交叉 相关 方法 来 确定 DA 白矮星 的 APP 速度 ， 再 减 去 白矮星 的 引力 红 移 就 得 到 了 白矮星 的 视 向 
速度 。 测 试 发 现 对 于 Teff 高 于 10000K 的 DA 白矮星 的 低 分 辨 率 光 谱 (R 2000) ， 在 信 噪 比 大 于 20 的 情况 
下 ， o E E 0 
视 向 速度 的 平均 值 接近 


关键 字 白 矮星 ， 视 向 速度 ， 模 板 匹 配 
中 图 分 类 号 : P141.5 


1 引言 


白矮星 是 银河 系 绝 大 部 分 恒星 的 最 终归 宿 ， 银 河 系 中 97% 的 恒星 最 终 都 要 演化 成 白矮星 。 白 矮星 由 于 
内 部 停止 了 核反应 ， 它 的 演化 非常 简单 ， 就 是 不 断 的 冷却 。 温 度 越 来 越 低 ， 光 度 越 来 越 暗 ， 最 后 变 成 一 颗 
黑 矮 星 。 根 据 光 谱 型 ， 白 矮星 可 以 分 为 DA、DB、DC、D0、DZ、DQ 等 类 型 ， 其 中 DA 白矮星 是 最 多 的 ， 大 
约 占 总 数 的 75%。 由 于 白矮星 光度 低 ， 观 测 相对 困难 ， 特 别 是 获得 白矮星 的 光谱 尤其 困难 。 在 Sloan 数字 
巡天 包 之 前 ， 最 大 的 光谱 证 认 的 白矮星 星 表 是 MaCcook 等 人 外 在 1999 年 发 布 的 ， 总 数 仅 有 2000 颗 。 在 他 
peg WINS EGE YS APRIL, ES, WA, MATAR, MPENDA RERU, WAER AA 
Ai 星 的 大 气 参数 。Sloan MAM K 1 tHE 4 AREWA ERK S 10 倍 以 上 。2006 4, Eisenstein 等 
>= A“ AE SDSS DR4 中 证 认 出 9000 多 颗 白矮星 ，Kleinman 等 人 外 在 SDSS DR7 中 证 认 出 了 20000 多 颗 白 矮星 ， 
= Kepler 等 人 中 2015 年 在 SDSS DR10 中 也 发 现 了 9000 多 个 新 的 白矮星 。 目 前 光谱 识别 的 白矮星 样本 已 经 超 
过 了 30000 多 个 ， 这 对 于 研究 白矮星 的 统计 特性 提供 了 数据 基础 。 
© 运动 学 可 以 帮助 我 们 理解 白矮星 的 星 族 ， 因 此 白矮星 的 视 向 速度 测量 是 非常 重要 的 。 白 矮星 的 视 向 速 
= E RA APP 速度 ) 测量 相对 较 少 ， 主 要 原因 是 以 前 的 白矮星 光谱 的 信 噪 比 和 分 辩 率 较 低 。 最 开始 的 白 矮 
星 运动 学 研究 主要 依靠 白矮星 的 自行 ， 而 简单 的 把 视 向 速度 设置 为 零 。 首 次 比较 精确 的 测量 白矮星 视 向 
速度 的 工作 是 Pauli 等 人 ”外 ,但 是 白矮星 光谱 的 分 辩 率 比较 高 。 测量 低 分 辩 率 DA 白矮星 光谱 的 视 向 速度 
是 一 个 非常 困难 的 工作 。DA 白矮星 光谱 的 特点 就 是 在 光学 波段 主要 由 比较 宽 的 氢 线 和 连续 谱 构 成 ， 另 外 由 
于 stark 效应 ， 使 得 氢 线 的 轮廓 是 不 对 称 的 中， 这 使 得 白矮星 的 视 向 速度 测量 误差 较 大 。 为 了 得 到 低 分 辨 
率 白 矮星 光谱 的 视 向 速度 ， 本 文 介绍 了 一 种 基于 白矮星 参数 来 选取 模板 光谱 通过 相关 方法 来 获得 DA 白矮星 
的 视 向 速度 。 由 于 白矮星 除了 多 普 勒 效应 外 ， 还 有 引力 红 移 效应 ， 所 以 我 们 这 个 将 白矮星 的 总 体 效应 对 应 
为 APP 速度 。 
第 二 章 详细 介绍 了 我 们 的 方法 。 第 三 章 介绍 方法 对 于 模拟 数据 的 应 用 效果 。 第 四 章 介 绍 了 用 这 种 方法 
来 计算 一 些 实际 观测 光谱 的 应 用 效果 。 第 五 章 对 这 种 方法 进行 了 总 结 。 


ou 


2 基于 白矮星 参数 确定 模板 的 交叉 相关 方法 确定 DA 白矮星 的 APP 速度 
2.1 DA 白矮星 光谱 
DA 白矮星 光谱 非常 有 特点 ， 也 比较 好 辨认 ， 特 别 是 温度 比较 高 的 星 。 图 一 上 展示 了 一 颗 典 型 的 
白矮星 和 一 颗 主 序 星 的 光谱 。 绿 色 实 线 是 一 颗 主 序 星 的 光谱 ， 黑 色 实 线 是 一 颗 白 矮 星 的 光谱 。 白 矮星 的 四 
a 资金 项 目 : 国家 自然 科学 基金 (11703051); 973 M H (2014CB845704) 
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HB. Hy. and H6 在 下 面 用 黑色 的 实 线 标 出 。 可 以 看 出 DA 白矮星 光谱 的 特点 是 线 非常 少 ， 


仅 有 几 条 比较 强 的 氢 线 ， 而 且 线 特别 宽 。 这 两 条 光谱 具有 相同 的 分 辩 率 ， 但 是 就 同一 条 和气 线 相 比 ， 主 序 星 
的 氧 线 要 比 白矮星 的 氧 线 窗 很 多 。 这 导致 如 果 仅 依靠 一 条 氢 线 来 来 测量 视 向 速度 ， 由 于 宽 线 的 线 心 不 容易 


确定 ， 误 差 会 相对 较 大 。 所 以 ) 


] 氧 线 线 心 来 估计 白矮星 的 APP 速度 效果 并 不 理想 。 
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矮星 光谱 ， 其 中 


1 上 图 : 绿色 实 线 是 一 颗 主 序 星光 谱 ， 黑 色 实 线 是 一 颗 白 和 
Figure 1 Top: Green solid line represents the spectrum of a Main S 
one white dwarf star. Ha, HB. Hy, and HS are marked with short so 


图 1 下 展示 了 


100 200 


RARE FIRE. FR: 白矮星 光谱 的 四 条 氧 线 的 放大 显示 
equence star while the black solid line is the spectrum of 
id lines. Bottom: the zoomed Ha, HB, Hy. and Hô lines. 


矮星 的 四 条 氢 线 的 放大 图 。 从 上 到 下 分 别 是 ，Hau、HB 、HY 、and H6 。 这 四 条 线 都 表 


现 出 轮廓 的 非 对 称 性 ， 也 就 是 线 心 两 边 不 一 致 。 这 是 由 于 stack 效应 产生 的 。 
这 对 利用 线 心 来 估计 APP 速度 造成 了 一 定 的 困难 。 


HB. Hy, and HS 呈现 明显 的 轮 廊 不 对 称 。 


2.2 模板 匹配 和 交叉 相关 方法 计算 DA 白矮星 的 APP RE 
计算 APP 速度 的 方法 主要 有 两 种 ， 一 种 是 利用 多 条 谱 

FMR RH. A 

基于 模板 光谱 的 交叉 相关 方法 来 估计 视 向 速度 。 

我 们 利用 白矮星 的 理论 光谱 来 作为 模板 光谱 。 模 板 光 谱 是 


Ha 相对 好 一 些 ， 


线 的 线 心 来 估计 ， 另 外 一 种 就 是 利用 模板 光谱 基 


矮星 的 光谱 特点 ， 利 用 线 心 方法 来 估计 不 能 得 到 准确 的 速度 ， 所 以 我 们 利用 


Koester 计算 6 的。 整个 模型 光谱 库 的 有 效 温度 


范围 从 5000K 到 80000K, 表面 重力 加 速度 的 范围 从 7.0 到 9.5。 


为 了 得 到 准确 的 APP 速度 ， 首 先 要 从 模型 光谱 库 中 选 出 最 


适合 的 模板 光谱 。 这 里 说 的 适合 ， 就 是 选 出 的 模 


可 


板 光 谱 要 和 目标 光谱 上 共有 较 一 致 的 大 气 参数 ， 即 相近 的 有 效 温度 和 表面 重力 加 速度 。 


模板 匹配 过 程 如 下 : 首先 对 模型 光谱 和 观测 光谱 进行 


每 一 条 理论 光谱 和 观测 光谱 的 流量 比 ， 


在 没有 谱 线 的 位 置 选 出 8 个 连续 谱 窗口 (3862,3939,4035,4219,4490,4693,5025,5254] 。 在 这 8 个 窗口 计算 
然后 拟 合流 量 比 和 波长 的 5 阶 多 项 式 函 数 。 随 后 整 条 理论 光谱 和 这 


重新 采样 ， 间 隔 为 0.1A， 根 据 DA 光谱 的 特点 ， 


个 多 项 式 相 乘 并 计算 相 乘 多 项 式 后 的 理论 光谱 和 观测 光谱 的 卡 方差 。 最 小 卡 方差 对 应 的 理论 光谱 作为 我 们 


Br 


最 终 选取 的 模型 光谱 。 


选 出 模型 光谱 后 ， 我 们 用 交叉 相关 方法 来 估计 和 白矮星 的 APP 速度 。 流 程 如 下 : 


首先 设 定 最 小 速度 Vmin 最 大 速度 V_max 和 速度 增 量 V_delt。 然 后 从 最 小 速度 开始 ， 


目标 光谱 进行 相应 的 


移动 后 和 模板 光谱 计算 相关 系数 ， 然 后 增加 一 个 速度 增 量 继续 计算 ， 一 直到 最 大 速度 ， 这 样 就 建立 了 一 个 
相关 系数 和 速度 的 关系 ， 也 就 是 相关 函数 。 接 下 来 寻找 相关 函数 的 最 大 值 ， 最 大 值 对 应 的 速度 就 是 我 们 最 


终 计算 得 到 的 APP 速度 。 
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2.3 方法 测试 


为 了 测试 方法 的 可 靠 性 ， 我 们 生成 了 一 
噪 比 的 光谱 。 


噪音 ， 来 形成 一 定 信 


测试 的 结果 。 


RV (80 km/s) 


系列 的 模拟 光谱 。 
然后 再 假设 一 定 的 APP 速度 ， 进 行 相应 的 移动 ， 作 为 一 个 模拟 光谱 。 然 
后 我 们 利用 2.2 的 方法 对 这 些 模拟 光谱 进行 了 APP 速度 的 测量 ， 图 2 展示 了 对 信 噪 比 为 20 的 模拟 光谱 进行 
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模拟 光谱 的 生成 是 在 理 六 


从 光谱 的 基础 上 增加 
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ZN 有 效 温 


图 1 测试 样本 的 误差 和 温度 的 关系 。 


Figure 1 Relation between the error and the e 


Ta (K) Ta (K) 


ective temperature 


图 2 中 ， 从 左上 角 到 右 下 角 展 示 了 APP 速度 为 80km/s ~ -80km/s 的 模拟 光谱 的 测试 结果 。 横 坐标 表 


度 ， 纵 坐标 表示 速度 的 误差 。 每 一 幅 


图 表示 一 条 误差 随 温度 变化 的 曲线 。 总 的 来 看 ， 误 差 在 低温 


区 相对 较 大 。 从 5000K 到 9000K， 误 差 呈现 下 降 趋势 。 在 5000K 的 时 候 误 差 最 大 。 而 且 总 体 误差 和 速度 


KAA 
于 10km/s. 


关 。 速 度 绝对 值 越 大 ， 误 差 相 对 越 大 。 对 于 有 效 温 度 高 了 


3 ”对 观测 样本 进行 APP 速度 测量 


我 们 利用 第 二 章 介 2 


包含 2 万 多 颗 白 矮星 ， 并 且 提 供 
的 DA 白矮星 大 约 2000 4. A 
相对 于 APP 速度 的 函数 ， 虚 线 指 的 是 最 大 相关 系数 对 应 的 速度 值 。 


星 的 APP 速度 。 


的 方法 对 SDSS 的 DR7 呈 中 的 DA 白矮星 样本 进行 了 APP 速度 测量 。 ) 
1 效 温度 ， 表 面 重力 加 速度 ， 质 量 等 


了 白矮星 的 


3 展示 了 对 一 


于 10000K 的 光谱 ， 白 矮星 的 速度 的 误差 小 


这 个 样本 共 

等 信息 。 其 中 信 品 比 大 于 20 
颗 DA 白矮星 进行 APP 速度 测量 的 例子 。 实 线 表示 相关 系数 
这 个 速度 就 是 我 们 计算 得 到 的 这 颗 白 矮 
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RV_APP=—19.2 km/s 


Corr 
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图 3: 计算 一 颗 白 矮星 APP 速度 的 相关 函数 ， 横 坐标 是 速度 序列 ， 纵 坐标 是 相关 系数 。 虚 线 指 的 是 最 大 相 

关系 数 对 应 的 速度 的 值 。 

Figure 3: The correlation function of APP velocity measurement for a DA white dwarf. Dashed line indicates 

the position of max correlation coefficient. 
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图 4: 观测 光谱 和 理论 模型 光谱 拟 合 的 一 个 例子 。 黑 色 实 线 为 观测 的 光谱 ， 红 色 需 线 为 相应 的 匹配 最 佳 理 
论 的 光谱 。 
Figure4: An example for fitting between the observed spectrum and the theoretical spectrum. Black solid 
line represents the observed spectrum while the red dotted line is the theoretical spectrum. 
3.1 白矮星 的 引力 红 移 计算 

在 计算 APP 速度 之 前 ， 首 先 我 们 计算 一 下 白矮星 的 引力 红 移 。 公 式 如 下 : 


(1) 
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MM 表示 白矮星 的 质量 ，10gg9 表示 白矮星 的 表面 重力 加 速度 ，R 表 示 白 矮星 的 半径 ，logg 鲍 康 示 太阳 的 表面 重力 
加 速度 ， 这 里 我 们 取 4.44。 这 个 公式 来 自 于 万 有 引力 公式 的 变换 。 


Vg=0.635 (2) 
人 表示 和 白矮星 的 质量 ，R 表示 白矮星 的 半径 。vVg 表示 白矮星 的 引力 红 移 。 这 个 公式 来 自 于 Silvestri GAM 
的 工作 


首先 根据 白矮星 的 质量 来 根据 公式 〈1) 计算 白矮星 的 半径 ， 然 后 根据 公式 (2) 来 估计 和 白矮星 的 引力 红 移 。。 
图 5 的 上 图 显示 了 这 个 样本 的 引力 红 移 和 视 向 速度 的 分 布 ， 可 以 看 出 引力 红 移 主要 分 布 在 10-80km/s。 大 部 分 
集中 在 30km/s 附近 。 这 也 和 和 白矮星 的 质量 主要 分 布 在 0. 6 太阳 质量 一 致 。 
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Dis (pc) 
图 5 上 图 : 整个 观测 样本 的 白矮星 的 引力 红 移 和 视 向 速度 的 关系 下 图 : 观测 样本 中 白矮星 的 视 向 速度 和 日 心 距离 的 关系 ， 绿 色 虚 线 表 示 视 向 速 
ERKO 的 线 ， 紫 色 和 矩形 框 表示 某 一 日 心 距离 范围 内 的 视 向 速度 的 平均 值 及 其 弥散 
Figure 5 Top: relation between the gravitational redshift and the radial velocity; Bottom: relation between the radial 
velocity and the distance. Rectangles represent the average radial velocity with a distance bin. 


3.2 白矮星 的 视 向 速度 测量 
对 于 这 2000 个 白矮星 样本 ， 我 们 测量 了 白矮星 的 APP 速度 。APP 速度 和 引力 红 移 之 差 就 是 白矮星 的 
视 向 速度 。 图 4 展示 了 确定 白矮星 模板 光谱 的 例子 。 红 色 虚 线 为 相对 于 观测 光谱 (黑色 实 线 ) 符合 最 好 的 


模板 光谱 ， 这 个 模板 光谱 用 来 做 相关 来 计算 
心 距离 的 关系 。 我 们 统计 了 日 


矮星 的 


lB 100pc. FAFA 


是 合理 的 。 
4 结论 


基于 白矮星 参数 来 选取 模板 进行 
白矮星 ， 能 够 得 到 较 准确 的 视 向 速度 ， 


匡 表 示 每 一 个 区 间 内 视 向 
离 在 100pc 和 1000pc 之 间 ， 视 向 速度 平均 值 在 0 PAE. IK Hil 


速度 的 


APP 速度 。 图 5 的 下 图 展示 了 
心 距离 在 1000pc 范围 内 的 视 向 速度 的 平 


匹配 来 测量 DA 


速 


矮星 的 视 向 速度 和 日 
均值 和 弥散 。 共 分 为 了 10 个 
平均 值 ， 误 差 棒 用 弥散 来 表示 。 可 以 看 出 ， 在 
FE 实 了 以 前 的 对 白矮星 的 视 向 速度 假设 为 


区 间 ， 
日 心 距 


0 


矮星 的 视 向 速度 ， 对 于 有 效 温 度 高 于 10000K 的 DA 
精度 在 10km/s 以 内 。 通 过 对 2000 多 颗 高 信 品 比 的 白矮星 的 视 向 


度 分 析 ， 我 们 发 现 视 向 速度 平均 值 整体 趋向 于 0， 这 表明 以 前 的 研究 在 计算 白矮星 的 运动 学 信息 中 假设 


矮星 的 速度 为 0 是 比较 合理 的 。 
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Radial Velocity Measurement of the DA White Dwarfs 
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AbstractWhite Dwarfs (WDs) are the final stage for the evolution of the majority of low- and medium-mass 
stars with initial mass < 8 Mo. Approximately 75% of all observed WDs are DA type. The spectra of DA white 
dwarf are dominated with Balmer Hydrogen lines. It is very difficult to measure the APP radial velocity of WDs by 
the line center because of the wide and asymmetry spectral profile. The kinematics of WDs is very important to 
distinguish populations of WDs. This paper present a method for APP radial velocity measurement using cross- 
correlation between the observed spectrum and the template spectrum that has the similar parameters with observed 
spectrum, which has advantage than linecenter method and the template spectra are from theoretical spectra by 
Koster et al. The final radial velocity is derived by the difference of APP radial velocity and the gravitational 
redshift. After test with simulated spectrum, the accuracy will be within 10 km/s for DA white dwarfs with Torr 
higher than 10000 K and SNR larger than 20. In addition, the radial velocities of observed DA WDs from SDSS 
DR7 are measured with the above method. The average radial velocities of these close DA white dwarfs (distance < 
1000 pc) are approaching zero . 


Keywordwhite dwarf; radial velocity;template match; 


